Low frequency noise controls on-off intermittency of bifurcating systems 
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A bifurcating system subject to multiplicative noise can display on-off intermittency. Using a 
canonical example, we investigate the extreme sensitivity of the intermittent behavior to the nature 
of the noise. Through a perturbative expansion and numerical studies of the probability density 
function of the unstable mode, we show that intermittency is controlled by the ratio between the 
departure from onset and the value of the noise spectrum at zero frequency. Reducing the noise 
spectrum at zero frequency shrinks the intermittency regime drastically. This effect also modifies the 
distribution of the duration that the system spends in the off phase. Mechanisms and applications 
to more complex bifurcating systems are discussed. 

PACS numbers: 05.45.-a, 91.25.Cw 



Among the possible behaviors of a chaotic system is 
intermittent behavior. The system remains for long du- 
rations in some regular state (say a laminar state or off- 
phase) and at unpredictable instants begins to explore 
other states (say on-phase) before returning to the lami- 
nar state. 

A simple deterministic model for intermittency was 
proposed by Pomeau and Manneville pj: a limit cycle 
is weakly unstable but from time to time a reinjection 
mechanism forces the system to return close to this limit 
cycle. A few years later, a new type of intermittency 
was discovered in coupled dynamical systems and also 
identified in a system of reaction-diffusion equations Q . 
Both systems can be approximately described by the evo- 
lution of a weakly linearly unstable mode with a noisy 
control parameter. This type of intermittency was given 
the name "on-off intermittency" by Piatt, Spiegel and 
Tresser Q who pointed out its genericity when an un- 
stable system is coupled to a system that evolves in an 
unpredictable manner. Experimentally, on-off intermit- 
tency has been identified in various systems including 
electronic devices, electrohydrodynamic convection in ne- 
matics, gas discharge plasmas and spin-wave instabilities 

a. ^ 

It is surprising that, despite the genericity of the on- 
off intermittency mechanism, this effect has not been re- 
ported more often. One might expect that any careful 
experimental investigation of an instability should reveal 
on-off intermittency when the system is close to the on- 
set of instability, and is hence sensitive to unavoidable 
experimental noise in the control parameters. This re- 
mark is the main motivation for the present work. We 
show that the amplitude of the noise is not the relevant 
control parameter of on-off intermittency. 

Through an analytical study of a simple stochastic sys- 
tem, we identify the parameter that drives the intermit- 
tent behavior and compare our prediction to numerical 
simulations. We then test our prediction for a chaotic 
rather than a stochastic system. Then we discuss the 



sensitivity of the statistics of the laminar phase duration 
to the parameter that controls the on-off intermittency. 
Finally, we present applications of this result to complex 
systems. 

One of the simplest systems that can exhibit on-off 
intermittency is 

X = (a + ((t))X-X\ (1) 

where £ is a random process with zero mean 0. In the 
deterministic regime (no-noise), the variable X under- 
goes a pitchfork supercritical bifurcation for a = 0. The 
attractor, X = 0, is stable for negative a and is unstable 
for positive a: X tends in the long time limit to one of its 
two stable attractors ±y/a. In the stochastic regime, the 
noise £ acts as a modulation of the forcing parameter. 
Note that if the initial condition verifies X(t = 0) > 
then X(t) > for all time. Henceforth we consider only 
positive initial conditions for X without any loss of gen- 
erality. 

For stationary Gaussian white noise, the probability 
density function (PDF) of X is derived by solving the as- 
sociated Fokker-Planck equation |(j . We define the noise 
intensity by (((t)((t')) s — D5(t — t') where () s is the av- 
erage over realizations of the noise. Equation Q is then 
understood in the sense of Stratonovich. If a < 0, X 
tends to zero and the stationary PDF is P(X) — S(X). 
For positive a, one obtains 



P(X) = CX^~ x e 



(2) 



where C is a normalization constant. For < 2a /D < 1, 
this PDF diverges at the origin X — 0. As pointed out 
in 0, this divergence is associated with the intermittent 
behavior of A", as A remains for long durations arbitrarily 
close to the unstable fixed point X = 0. When 2a/ D is 
large, intermittency disappears and X fluctuates around 
its deterministic value, y/a. 

However, for colored noise, a more complex situation 
is expected. It is tempting to assume that the noise am- 
plitude v/ (C(t) 2 ) s controls the on-off regime. This is not 
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FIG. 1: Temporal traces of the dynamical variable X(t), 
solution of JTJ, with a = 1.25 x 10~ 3 , a 2 = 0.005 and the 
noise defined by ©: (a) r) = fi = 0.25 i.e., a/S = 0.3927; (b) 
77 = O = 2.5 i.e., a/5" = 3.9270. 



the case. We plot in fig. [I] the solution of equation for 
two different colored noises with the same value of (C 2 )s- 
One of the solutions is intermittent but not the other. 
Hence intermittency is controlled by another parameter 
of the system. 

In order to identify this parameter, we derive an ap- 
proximate expression of the stationary PDF of X using 
the cumulant expansion introduced by van Kampen and 
shown to be valid for small values of ar c where a is the 
noise amplitude and t c its correlation time Q. In the 
case under study, two parameters appear in the expan- 
sion that are related to the autocorrelation function: 



8 = 



(C(0)C(r)).dr , M= / (C(O)C(r)) 



Z- 2aT dT . 



(3) 

For a Gaussian white-noise, the expansion is exact and 
leads to the stationary PDF given by equation (J2J. In 
the following we consider the generic case where S and 
M — S are non-zero. The stationary PDF is 



P{X) = CXs-^l 



(M-S)X 2 (1+ 
S a 



2 S (M-S 



(4) 



where C is a normalization constant. The behavior of 

the PDF for X close to zero is proportional to 

It diverges for X = so that on-off intermittency occurs 

if 



0<-<l. 



(5) 



This is consistent with fig^ since a/S — 0.3927 for the 
intermittent signal and a/S — 3.927 for the other one. 
The Wiener-Khintchin theorem states that the integral 
of the correlation function, 2S, is equal to the spectrum 
of the noise at zero frequency. Thus, another interpre- 
tation of the criterion J5J is that on-off intermittency is 
present when the departure from the deterministic on- 
set is smaller than half the value of the noise spectrum 



at zero frequency. This interpretation also holds when 
the noise is white and Gaussian since S = D/2. For 
white-noise, the spectrum has the same value, D, for all 
the frequencies. Consequently, the analytical result 
does not identify which part of the spectrum controls the 
intermittency. 

In order to check the validity of this result, we solve 
equation Q numerically, using a stationary Gaussian 
correlated noise ((t) with autocorrelation function Q : 



(C(*)C(* + T ))s = a 2 (cos(2tt^t) + g sin(27rfi|r|)) 

exp (— 27t?7|t|) . 



(6) 



The noise variance is a 2 and its correlation time t c = 
(2ttt])~ 1 . This provides S — a 2 r]/ [Tr(r] 2 + ft 2 )] and 
M = a 2 (r) + o/(2tt))/ [tt ((77 + a/71-) 2 + ft 2 )] . Therefore 
by changing rj and ft we can tune a/S and ar c inde- 
pendently. Gaussian white noise is recovered in the limit 
7/ — » 00 with a 2 /i] = D. Figure ^ presents time series of 
X and fig. [3 the corresponding PDFs. We also compare 
the predicted criterion for appearance of intermittency 
(jSJl with the numerical results. To wit, we draw in fig- 
ure |21 a phase diagram in the (S, a)-plane using noises 
with different S and M. We calculate X mp the most 
probable value of X. For a > 0, the solution X = is 



unstable. The system is intermittent if X. n 



0, and 



non-intermittent if X mp 7^ 0. In these figures and for all 
tested parameter values for which ar c is small, there is a 
very good agreement between the numerical results and 
the predictions (|4I5() . 




FIG. 2: PDF of the solutions of JTJ with the noise ©. The 
symbols (x) and (+) correspond, respectively, to the param- 
eters used in figures la and lb. The full lines are the corre- 
sponding predictions given by IpT). Same figure in inset using 
log-log scale. 

Up to now, we have only dealt with fluctuating param- 
eters that are random processes. It is tempting to test the 
prediction (1 11511 with a deterministic but chaotic fluctu- 
ating parameter. Thus, we study equation Q when £ is 
obtained from the chaotic solution of the Lorenz system 
131, i.e., we solve 



U 



-a(U-Y),Y = rU-Y-UZ,Z =UY-bZ, (7) 
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FIG. 3: Behavior of the solution of Q for noises defined in 
© and associated with various values of S and M. (A) : 
intermittent behavior, (o) : non-intermittent behavior. The 
full line is the transition curve predicted by JSJ. 



and define £ by 



liU n + (1 - n)U n 
! (( M [/ n + (l- M )tf„)2)l/2 



(8) 



U„ = 



u-(u) 



and U n — 



u-(U) 



Averages are 



now understood as long time averages and y/ < £ 2 > = a. 
The parameter fi is tuned between zero and one to change 
the amplitude of the noise spectrum at zero frequency. 
Since U is the derivative of U, the value of its spectrum 
at low frequencies is smaller than that of U. Increasing 
H increases the weight of U and thus decreases the noise 
spectrum at low frequencies (accordingly the value of S). 

The equations Q 0) are then solved numerically for 
r = 25, a = 10 and b — 8/3. The solution of equation Q) 
is chaotic and we plot examples of time series of £ and 
X in fig. 0] As expected on-off intermittency disappears 
when /i increases, i.e., S decreases. For the intermittent 
signal we have a/S ~ 0.332, but a/S ~ 5.64 for the non- 
intermittent one. The numerical estimates of the PDFs of 
X are compared to the fit given by (0J (in fig. |3J). Again, 
for small values of the noise amplitude, the agreement 
between the prediction and the numerical results is very 
good. 

Interesting results can also be obtained for the dura- 
tion of the laminar phases r. We define a laminar phase 
as follows: X(t ) = e, X(t) < e for t < t < t + r 
and X(to + r) = e, e being an arbitrary threshold below 
which the system is considered to be in the laminar state. 
Close to the onset of on-off intermittency, the probabil- 
ity P(t) of the duration t of the laminar phase satisfies 
P(t) oc t~ 3 / 2 [nj]. For large values of r, a cut-off in the 
power law appears at finite departure from onset . We 
have checked numerically that for a colored noise, the 
PDF is indeed proportional to r~ 3 / 2 with a cut-off for 
high r. The position of the cut-off increases when S in- 
creases (a being constant). This is consistent again with 
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FIG. 4: Bottom panels: solutions of with ( obtained from 
the Lorenz system ©. The chaotic functions £ are displayed 
in the top panels. The parameters are a — 0.01, ^/ < C, 2 > = 
0.2 and fj, = (left panels), [i = 0.8 (right panels). Note the 
difference in the time scales t between the top and bottom 
panels. 
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FIG. 5: PDF of the solutions of Q with ( obtained from 
the Lorenz system JSJ. The parameters are the same as in 
figure ^] with (x) /i = 0, (+) /i = 0.8. The inset is the same 
figure with log-log scales. The full lines are the corresponding 
predictions (gj. 



our interpretation of the role of the noise spectrum at 
zero frequency: the smaller a/S is, the longer the system 
remains in the laminar state and the more intermittent 
the signal appears. 

Our interpretation of the phenomenon is as follows. 
On-off intermittency occurs because of a competition be- 
tween the noise and a systematic drift driven by the dis- 
tance from onset. More precisely, as pointed out in 0, 
when X is close to the unstable manifold X = 0, the 
evolution of Y = log A is given by Y = a + ((t). For 
positive a, Y is positive on average but events in which 
Y remains smaller than its initial value are possible pro- 
vided / = / ((t)dt/T remains smaller than —a for a long 
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duration. In the long time limit, the main contribution 
to the integral I is due to the zero frequency component 
of the noise. If this component is reduced, then occur- 
rences of the inequality I < —a are less probable and 
intermittency tends to disappear. Note that our ana- 
lytical calculations are based on perturbative expansions 
and are valid for small values of the product of the noise 
amplitude with its correlation time. However, even for a 
finite amplitude of noise, we have verified that when the 
low frequencies arc filtered out, intermittency disappears. 

We expect that the role of the zero frequency compo- 
nent of the noise is generic and also pertinent for systems 
more complex than the one presented here. For instance, 
Liicke and Schanck studied a system similar to eq.(|T|l 
with inertia taken into account. Through a perturba- 
tive expansion close to the deterministic solution, they 
calculated a noise-induced postponement of the onset of 
instability and a modification in the amplitude of the un- 
stable mode. As they pointed out later, their expansion 
is not correct when the noise spectrum does not vanish 
at zero frequency Recent calculations for the same 
system have shown that for a small departure from onset 
and for an Ornstein-Uhlenbeck or a white noise, the PDF 
of the unstable mode diverges close to zero hj. In the 
light of our study, both the divergence of the PDF and 
the failure of the perturbative expansion are related to 
the same physical effect: on-off intermittency when the 
noise spectrum at zero frequency is non-zero. 

Our analysis explains why many experimental investi- 
gations on the effect of a multiplicative noise on an in- 
stability do not display on-off intermittency. If the noise 
is high-pass filtered, as often required for experimental 



reasons, then the regime of intermittent behavior disap- 
pears. This is the case for instance in a ferrofluidic 
layer undergoes the Rosensweig instability and peaks ap- 
pear at the surface. The layer is then subject to a mul- 
tiplicative noise through random vertical shaking. Close 
to the deterministic onset, the unstable mode submitted 
to a colored noise does not display intermittency. 

In dynamo theory, the magnetic field is forced by the 
flow of an electrically conducting fluid. The velocity of 
the flow appears as a multiplicative term in the equation 
for the magnetic field. If the flow topology is complex 
enough and the velocity is large, a magnetic field is gener- 
ated by dynamo instability. The flow is in general turbu- 
lent at dynamo onset so that the velocity fluctuates. We 
infer that the intermittent behaviors as seen numerically 
by Sweet et al. 0] are related to the presence of very 
low frequencies in the spectrum of the velocity field. The 
same features also occur in simple models of dynamos 
subject to white noise 0|. On the contrary, experimen- 
tal realizations of dynamos driven by constrained flows 
did not display intermittency 01- These flows, though 
turbulent, are probably too constrained to display ve- 
locity fluctuations with large enough amplitude of the 
spectrum at low frequencies. 

This work could be generalized to study a parametric 
instability with a time-dependent forcing. For an har- 
monic forcing subject to frequency or amplitude noise, 
intermittent behaviors have been reported |18| . In these 
cases, the relevent component of the noise that controls 
on-off intermittency still remains to be identified. 

This work has greatly benefited from discussions with 
N. Leprovost, N. Mordant, J. Farago, S. G. Llewellyn 
Smith, S. Fauve and P. Marcq. 
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